Pattern formation upon femtosecond laser ablation of transparent dielectrics 
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Costache et al. have reported recently a new type of periodic patterns generated at femtosecond 
laser ablation of transparent dielectrics (Appl. Surf. Sci. 186, 352(2002)). They show features 
known from other pattern forming systems far from equilibrium, like point and line defects or grain 
boundaries, and cannot be explained by the classical theory. The present work is an attempt to 
investigate these pattern by means of a generalized Kuramoto-Shivashinsky equation derived from 
the Bradley et al. and Cuerno et al. model for ripple formation at ion beam sputtering of surfaces. 



Laser induced periodic surface structures (LIPSS) have 
been observed for almost 40 years [l[ on targets made 
of intrinsic and extrinsic semiconductors, metals and di- 
electrics, using cw to subpicosecond laser sources with 
wave lengths varying from the ultraviolet up to the in- 
frared domain (0 and references therein, 0,0,0). They 
have therefore been considered to be a universal phe- 
nomenon that can occur on any material that absorbs 
radiation, regardless of its dielectric constant . 

The explanation accepted today on a large 
scale has been delivered by the group of 
ISipe. Young. Preston, and van Driel in Il983l by 
taking into account the details of the interaction of an 
electromagnetic wave with the microscopically rough 
selvedge of a surface. Thus, the periodic surface struc- 
tures are the result of the interference between the 
incoming laser light, the light refracted by the bulk and 
the light scattered by the rough selvedge (or interface). 
The aspect of the surface is usually a combination 
of several basic patterns. For p-polarized light, these 
are the s-type patterns (Eq. Q)) perpendicular to the 
electric field E and the c-type patterns parallel to E; for 
s-polarized light, there are c-type patterns perpendicular 
to E and other patterns with no simple dependence 
on the angle of incidence. Excellent agreement with 
experimental work that followed 0] has been found. 
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A denotes here the wave length of the laser and A the 
periodicity of the patterns. 

Recently, ICostache et alJ have reported upon a new 
type of periodic patterns generated at femtosecond laser 
ablation of transparent dielectrics 6]. These new struc- 
tures exhibit fundamental differences when compared to 
the "classical" periodic patterns, such as spacing almost 
twice as small, independent on the wave length and on 
the angle of incidence of the laser beam, but correlated 
with the local incident intensity. Features like point and 
line defects or grain boundaries can also not be described 



The qualitative and quantitative similarities to the 
patterns observed in ion beam sputtering 0, as 
well as the wide accepta nce of their theoret ical de- 
scription as proposed by iBradlev and Harper] and 
ICuerno and Barabasil , have determined us to use this 
framework for studying the new structures. 

The experiments 0, [HI E| have been carried out 
under high vacuum (< 10 _T mbar) on freshly cleaved sin- 
gle crystal slides of BaF2 and CaF2. The laser system 
generated a pulse of 120 fs with a central wave length 
of 800 nm and intensities of 1 - 12 x 10 12 W/cm 2 . Ad- 
ditionally, the frequency of the laser radiation could be 
doubled and the angle of incidence of the laser beam on 
the target could be varied. 

The main feature which distinguishes these new struc- 
tures from the classical ones is the spacing. In the clas- 
sical model the patterns are uniformly distributed across 
the ablation area, with spacings given by Eq. (QJ. In 
this case, however, no dependency of the ripples period- 
icity on the laser wave length or on the angle of incidence 
could be observed. Moreover, the periodicity is obviously 
dependent on the local intensity of the laser beam, as 
shown in Fig. |l(a)| The spacing of the ripples increases 
here towards the center of the ablation spot, that is, in- 
creases with the local intensity. While Eq. Q allows a 
minimum spacing of w 468 nm for the situation shown 
in Fig. |l(a)[ the spacing at the boundary of the ablation 
spot is approximatively 250 nm. 



The model prese nted bv IBra dlev and H arperl and 
further developed bv lCuerno and Barabasil \^ is based on 
Sigmund's theory of sputtering [l^. Here, an ion striking 
a solid will first travel a certain distance a, called average 
depth of energy deposition, and then lose its energy in a 
cascade of random atomic collisions. For an ion travel- 
ling along the z-axis, the resulting profile of the average 
energy deposition in the material (z < h(x, y)) follows in 




FIG. 1: Orientation of the periodic ripple structure with respect to the laser beam polarization, (a) Ablated spot obtained 
after 9200 shots at 0.8 x 10 13 W/cm 2 ; (c) 5000 shots at 1.2 x 10 13 W/cm 2 ; (a), (b) and (c) 800 nm, 45° incidence. 
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Fd{t) denotes here the energy deposition per unit vol- 
ume, e is the total energy deposited by the ion and a 
and j3 are the widths of the distribution parallel and per- 
pendicular to the plane of incidence, respectively. The 
erosion velocity at a certain point on the surface is then 
proportional to the total energy deposited there by the 
ions within the range 1Z of the distribution @ 
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where \I/(r) is the local flux of ions and A is a material 
constant. 

Eq. (|2J) was written in the reference frame of the in- 
coming ion, having the origin at the point of impact O 
and the z-axis pointing away from the surface, along the 
velocity v of the ion. In Eq. Q we went over to the local 
reference frame, with the origin also at O but with the z- 
axis pointing along the surface normal. And finally, there 
is the reference frame of the laboratory, our destination. 

To evaluate the erosion velocity we make the trans- 
form x — > a( x , y — > a( y and expand the integrand in 
powers of a/Rx and a/Ry- Rx and Ry are the curva- 
ture radii at O and are used to describe the integration 
region 1Z 
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with 1/Rx = -d 2 h/dx 2 and 1/R Y = -d 2 h/dy 2 . After 
evaluating the resulting Gaussian integrals we make the 
transition to the laboratory coordinate frame (x, y, h) || 
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and obtain thus for a = /? ICnerno and Barabasi ls equa- 
tion of motion y| 
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In Eq. (0 (j) denotes the angle between the incident beam 
and the surface normal at O. It can be expressed as a 
function of the angle of incidence and the local gradient 
Vh and can be expanded in powers of the latter. 

Eq. © has the form of the Kuramoto-Shivashinsky 
(KS) equation |l4j. i>o, /i, v x , v y , \ x and \ y are func- 
tions of the angle of incidence 0. The term —KV 2 (V 2 h) 
accounts for the surface self-diffusion and rj(x, y, t) is a 
Gaussian white noise, accounting for the stochastic ar- 
rival of the ions on the surface. If the surface diffusion is 
thermally activated, the coefficient K is given by [jj 



K = Ds^v '/n 2 kBT, Ds 
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where Ds [l5( is the surface self-diffusivity, Q a is the 
activation energy, 7 is the surface free energy per unit 
area and v is the areal density of diffusing atoms. 
In the linear regime, periodical stripes with 
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have the highest growth rate. Consequently, for v = v x < 
the surface will be dominated by ripples along the y- 
axis (k || x) and for v — v y < by ripples parallel to the 
x-axis (k || y). 

In the case of the laser material interaction, we assume 
that the photons loose their energy in a similar stochastic 
process. The widths a and (3 of the Gauss distribution 
(|2j are equal because the target is isotropic. The inci- 
dent energy flux e\P(r) is now given by aj\S2n\, where 
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a has the meaning of an average penetration depth, 7 is 
the absorption coefficient of the material and §2 is the 
Poynting vector of the refracted beam. The product 
represents the fraction of the transmitted energy that will 
be absorbed, that is, in an ionic picture, scattered by the 
material. Using the Fresnel equations, the incident flux 
can be written in terms of the laser intensity Io = \Si\ 
and of the transmission coefficient T 
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In the next step we insert the energy flux © and the 
Gaussian distribution @ into Eq. © to compute the 
local erosion rate. To evaluate the integral we extend 
the Taylor expansion Q to its complete form 
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where 1/Rxy = —d 2 h/dxdy. 

We call the last term in Eq. (IT0|) a "rotational" co rrec- 
tion. In their original model, iBradlev and Harper! have 
neglected it and could therefore only describe patterns 
growing along the x- or the y-axis. With this term 
the stability analysis can still be reduced to the case of 
IBradlev and Harper! but in a reference frame rotated by 
an angle B to the original one. Following the notation 
style of Eq. © and denoting by v xy the coefficient of h xy , 
the angle is determined by 



v xy cos 20 = [y x — v y ) sin 20. 
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We perform the transition ® to the coordinate frame 
of the laboratory and scale the time and the space ac- 
cording to 



x = ax \ y — ay , t 
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to obtain 
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The tuples ijklm in Vijkim describe the product 
hxhyhxxhyyh™y these coefficients belong to. Except for 
5, all the coefficients in Eq. ([T3|) have the form {,} = 
aexp[— ^- cos 2 #]/{,} (w, with a = a/a. 

Fig. [21 shows a 3D plot of — v X: —v y and —v xy against 
the angle of incidence 6 and the polarization cp. We chose 
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FIG. 2: Plot of v x , v y and v xy against the angle of incidence 
and the polarization ip. a = 4, n = 1.47 



<T=4, B = 1.2 • 10 -4 and n=1.47, the refraction index of 
BaF2 at A ~ 750 nm. As shown in the previous sec- 
tion, the larger of — v x > and —v y >0 determines the 
prevailing pattern and v xy will rotate this pattern by an 
angle (Eq. jDJ). ' 

We see that there is actually only a small region be- 
tween # m i n and #max, where the ripple orientation changes 
with the light polarization. In our case # m i n was 29.64° 
and #max = 34.04°. For angles smaller than m i n the rip- 
ples should always be perpendicular to the surface com- 
ponent of the laser beam, where as for angles larger than 
#max they should be parallel. The width of this region is 
usually between 1 and 5 degrees, increasing with increas- 
ing n. Its position varies slightly with cr, but strongly 
with n. Smaller n's and larger cr's will move it towards 
angles of up to 50 — 60°. 

The present work is an attempt to investigate the oc- 
currence of newly discovered laser induced periodic sur- 
face structures^! by means of the Bradley et al. [3 and 
Cuerno et al. [2| model. Based on Sigmund's stochastical 
theory of sputtering this model has managed to ex- 
plain in an unified framework most of the dynamic and 
scaling behaviors observed experimentally at ion bom- 
bardment of surfaces. It eventually traces the origin of 
ripple formation down to the instability caused by the 
competition between surface roughening (so called neg- 
ative surface tension) and surface diffusion (the positive 
surface tension). 

The use of the Fresnel equations and the inclusion of 
the v xy (d 2 h/ dxdy) term to the Kuramoto-Sivashinsky 
equation, have allowed us to confirm the experimentally 
observed countinous change of the ripple orientation with 
the laser polarization. However, our model predicts that 
this should only happen for angles of incidence confined 
in a small range between # m i n and # m ax- In our simu- 
lations this range corresponded to 29 — 34°, but could 
move up to 53 — 55°. 

The role of the local laser intensity can also be ex- 
plained. According to Eq. Q, for local temperatures 
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FIG. 3: Surface topography for 0=31.8° and <p=0°, 30°, 60° 
and 90° respectively (upper- left to lower- right). 



^locai below Q a , K is a monotonically increasing func- 
tion of Ti oca i and thus of the laser intensity 2i OC ai- A 
first estimation yields Q a > 0.53 eV |22j , much larger 
than the melting point of BaF 2 (1280 °C « 0.13 eV). 
The periodicity 27r/|fc| of the ripples (Eq. ((HJ) is thus an 
increasing function of ii OC ab just like in Fig. |l(a)| 

At normal incidence, all patterns with the wave vector 
satisfying \k\ = ^/—v/2B can occur with equal probabil- 
ity. The symmetry break has therefore been looked for 
and traces of it have been found in the 3 rd order terms. 
However, further investigation is needed, as we have not 
been able to simulate any ripples yet. 

It is well known that the spacing of the structures de- 
scribed by the KS equation increases with time according 
to A ~ t 1 . The value of 7 varies with different experimen- 
tal conditions 0, ^3 ano ^ different numerical simu- 
lations 0,0,^1- Considering that the exposure times 
varied across the experiments, an unknown value of 7 
makes the direct quantitative comparison of the numeri- 
cal simulations with the experiment impossible. A com- 
prehensive study of the time behavior of the structures is 
therefore a further requirement for their understanding. 

iFrost et all have recently shown that the compound 
nature of a material can considerably influence the scal- 
ing behavior of the surface structures in the early stages 
of ion sputtering 16]. Considering the extreme condi- 
tions of the laser ablation, the loss of atoms is with 
high probability different for Ba and F, so that a spa- 
cially varying concentrat ion of adatoms on the surface 
would arise. iMavr et al.l have shown that surface diffu- 
sion driven by a concentration gradient of the adatoms is 



able to gen erate structure coarsening j^j. According to 
IFrost et all this type of coarsening would be one of the 
causes for the deviation from the KS predicted scaling 
behavior they have observed in the early stages of Ar + 
sputtering of InP surfaces. 

Further experiments with rotating targets or circular 
polarized light could also confirm the common nature of 
this new type of LIPSS and of the ion sputtering pro- 
cesses. Studying the topography of simultan eously ro- 
tated and Ar + ion sputtered InP surfaces, IFrost et al.l 
have been able to find structures of highly hexagonal 
symmetry 

I. Georgescu would like to thank Mario de Menech for 
helpful discussions. 
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